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We discuss the application of wavelet transforms to a critical interface model, which is known to 
provide a good description of Barkhausen noise in soft ferromagnets. The two-dimensional version 
of the model (one-dimensional interface) is considered, mainly in the adiabatic limit of very slow 
driving. On length scales shorter than a crossover length (which grows with the strength of surface 
tension), the effective interface roughness exponent C is ~ 1.20, close to the expected value for 
the universality class of the quenched Edwards- Wilkinson model. We find that the waiting times 
between avalanches are fully uncorrelated, as the wavelet transform of their autocorrelations scales 
as white noise. Similarly, detrended size-size correlations give a white-noise wavelet transform. 
Consideration of finite driving rates, still deep within the intermittent regime, shows the wavelet 
transform of correlations scaling as 1//^'^ for intermediate frequencies. This behavior is ascribed to 
intra-avalanche correlations. 
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I. INTRODUCTION 

In this paper we use wavelet concepts [l|, 0, S| to dis- 
cuss assorted properties of a single-interface model which 
has been used in the description of Barkhausen "noise" 
(BN) [ijSSSBIil- BN is an intermittent phenomenon 
which reflects the dynamics of domain-wall motion in 
the central part of the hysteresis cycle in soft ferromag- 
nets (see Ref. [l3 for a review). By ramping an exter- 
nally applied magnetic fleld, one causes sudden turnings 
(avalanches) of groups of spins. The consequent changes 
in magnetic flux induce a time-dependent electromotive 
force V{t) on a coil wrapped around the sample. Analy- 
sis ofV{t), assisted by suitable theoretical modeling, pro- 
vides insight into both the domain structure itself and its 
dynamical behavior. It has been proposed that BN is an 
illustration of "self-organized criticality" 0, [ll|, [H, [l^ . 
in the sense that a broad distribution of scales (i.e. 
avalanche sizes) is found within a wide range of variation 
of the external parameter, namely the applied magnetic 
fleld, without any flne-tuning. The interface model stud- 
ied here 01 incorporates a self-regulating mechanism, in 
the form of a demagnetization factor. 

This way, real-space properties, e.g., interface rough- 
ness, reflect the divergence of the system's natural length 
scale, as it self-tunes its behavior to lie close to a second- 
order (interface depinning) transition. In this context, 
the application of wavelet transforms, which by construc- 
tion incorporate multiple length scales 0, 0, 13| , is natu- 
rally suggested. 

Also, when one considers the time series of intermittent 
events which characterizes BN, a broad range of varia- 
tion of V{t) is shown, in correspondence with the sim- 
ilarly wide distribution of avalanche sizes. Speciflcally 
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considering the model of Ref. H, it is known that the de- 
magnetizing term is responsible for the introduction of 
short-time negative (inter-avalanche) correlations (such 
correlations are observed in experiments as well) [4-, ^. 
Thus, a flnite time scale ("loading time") is introduced, 
which coexists alongside the broad distribution of V{t). 
The tool most frequently used in the analysis of BN time 
series is the Fourier power spectrum, i.e., the (cosine) 
Fourier transform of the time-time autocorrelation func- 
tion of the signal V{t) (lol. IT^. [l5| . BN power spectra ex- 
hibit distinct types of behavior along different frequency 
ranges, reflecting the fact that flnite "internal" times play 
relevant roles. For instance, the loading times referred to 
above are expected to influence the low-frequency end 
of the power spectrum, which pertains to inter-avalanche 
correlations, while the high-frequency tail relates to intra- 
avalanche ones. It has been stated that "understanding 
the power spectrum of the magnetization noise is a long 
standing problem" [l5| . 

Some existing appHcations of wavelet transforms to the 
analysis of V{t) jdJ, [13, [ll mainly aim at demonstrat- 
ing that the resulting spectra can successfully distinguish 
between BN originating from physically distinct mate- 
rials (e.g., samples under differing amounts of internal 
stress). Semi-empirical classiflcation schemes have been 
proposed (itI . Il8l |. Wavelet (Haar) transforms [l| have 
also been employed in conjunction with standard Fourier 
series, in order to produce hi gher -order power spectra 
of experimental data for V{t) [isl . [l9l . [2o| . Analysis of 
the corresponding results provides relevant evidence con- 
cerning correlations between events at different frequency 
scales. While in this work we shall deal only with flrst- 
order transforms, in Sec. [IV] below we shall comment 
on possible connections of our own flndings to those of 
Refs.Qlli,!!^. 

The paper is organized as follows. In Sec. [l] we recall 
pertinent aspects of the interface model used here, and of 
our calculational methods, as well as some basic features 
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of wavelet transforms. In Sec. IIIII we consider the scal- 
ing of interface roughness configurations. In Sec. IIVI we 
investigate properties extracted from time series, namely 
waiting-time and avalanche size correlations. Finally, in 
Sec. [Vl concluding remarks are made. 



II. MODEL AND WAVELET TRANSFORMS 
A. Single-interface model for BN 

We use the single-interface model introduced in Ref. Q 
for the description of BN. In line with experimental pro- 
cedure, the external field H acting on the sample is as- 
sumed to increase linearly in time, therefore its value is 
a measure of "time". Initially, we consider the adiabatic 
limit of a very slow driving rate, thus avalanches are con- 
sidered to be instantaneous (occurring at a fixed value 
of the external field). In this simplified version, a plot 
of V(t) against t consists of a series of spikes of varying 
sizes, placed at non-uniform intervals. Generalizations 
for a finite driving rate may be devised 1^, |2l|, [ISl; they 
are investigated in Subsection IIV Dl below. 

Simulations are performed on an L^; x Lj, x oo geometry, 
with the interface motion set along the infinite direction. 
Here we consider Ly — 1 (system dimensionality d = 
2, interface dimensionality d' = 1). Periodic boundary 
conditions are imposed at a: = 0, L . 

The interface (180-degree domain wall separating spins 
parallel to the external field from those antiparallel to it) 
is composed by L discrete elements whose x coordinates 
are Xi — i, i — 1,...,L, and whose (variable) heights 
above an arbitrary reference level are hi. The simulation 
starts with a fiat wall: hi = for all i. 

Each element i of the interface experiences a force given 
by: 



fi = u{xi,hi) + K [hi+i + hi 



2/1,: 



where 



He = H -r]M . 



(1) 



(2) 



The first term on the right-hand side of Eq. H]) represents 
quenched disorder, and is drawn from a Gaussian distri- 
bution of zero mean and width R; the intensity of surface 
tension is set by k, and the effective field is the sum of 
a time- varying, spatially uniform, external field H and a 
demagnetizing field which is taken to be proportional to 
M = (l/L) J2i=i the magnetization (per site) of the 
previously fiipped spins for a lattice of transverse width 
L. Here we mostly use R — 5.0, k = 1.0, 77 — 0.005, 
values for which fairly broad distributions of avalanche 
sizes are obtained jl, S, 0, Q • The exception is Sec. IIIH 
where (for reasons to be explained) , we allow the surface 
tension k to vary. 

The dynamics goes as follows. For fixed H, starting 
from zero, the sites are examined sequentially; at those 
for which fi > 0, hi is increased by one unit, with M 



being updated accordingly; the corresponding new value 
of u is drawn. The whole interface is swept as many times 
as necessary, until only sites with fi<0 are left, which 
marks the end of an avalanche. The external field is then 
increased until fi = for at least one site. This is the 
threshold of a new avalanche, which is triggered by the 
update of the site(s) with fi = 0, and so on. 

Because of the demagnetizing term, the effective field 
He at first rises linearly with applied field H, and then, 
upon further increase in H, saturates (apart from small 
fiuctuations) at a value rather close to the critical exter- 
nal field for the corresponding model without demagneti- 
zation 0, 01. 



B. Wavelets 

Wavelets are characterized by a scale parameter, a, and 
a translation parameter, b, such that the wavelet basis, 
{i^a-b{x)} can be entirely derived from a single function 
ij{x) through 



(3) 



The wavelet transform of a function f{x) is given by: 



W[/](a,6) = 



(4) 



Here we shall use the Daubechies wavelet family 
These are real functions (appropriate in the present case 
where the input signal is always a real number, whether 
it be an interface height or a voltage); in the discrete 
transform j2| implementation used here, the scales {a} 
are hierarchically distributed, i.e., aj = ag. We have 
experimented with the Daubechies wavelets of orders [l| 
4, 12, and 20, and found that, similary to Ref. 0, the 
quality of our results does not seem to depend on that. 
Therefore we have chosen the lowest order, Daub4, for 
our calculations. 

It must be noted that the Daubechies wavelet filter 
coefficients used here incorporate periodic boundary con- 
ditions In the applications to be discussed, for each 
case we shall comment on the specific consequences of 
this constraint. 

Furthermore, following Ref. 0, we have chosen to av- 
erage over the translation parameters 6, thereby arriving 
at a set of averaged wavelet coefficients to be denoted by 
14^[/](a). Among the several possible choices, we have 
found that averaging the squared coefficients tends to 
give smoother results than, e.g., using absolute values 0]. 
Thus, we define: 



W[f]ia)= {{W[fKa)f), 



1/2 



(5) 



where (•);, stands for average over the translation param- 
eters b. 
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III. REAL SPACE PROPERTIES: INTERFACE 
ROUGHNESS 

We begin by applying wavelet transforms to interface 
roughness data. The roughness W2 of a fluctuating in- 
terface with A'^ elements is the position-averaged square 
width of the interface height above an arbitrary reference 
level in,!!!: 



JV 



W2 = N' 



(6) 



where h is the average interface height. Self-afHnity prop- 
erties are expressed in the Hurst, or roughness exponent 

c MM- 



(7) 



where angular brackets stand for averages over the en- 
semble of allowed interface configurations, and [for the 
(1 -t- 1) dimensional systems which will be our main con- 
cern here] L is the profile length. 

Numerical evidence has been given [1] that, as regards 
interface configuration aspects, the model described here 
is in the quen ched Edwards- Wilkinson universality class. 
Thus [13, m, S [131 one expects C - 1.25 in d = 2 . 

We have simulated BN through the evolution in time 
of the adiabatic, d = 2 version of the model described 
above. Steady state, i.e., the stabilization oi He of Eq. ([2]) 
against external field H, occurs after some 200 events, 
for the range of parameters used here. In order to avoid 
start-up effects, here and in all subsequent sections we 
have used only steady-state data in our statistics. At 
the end of each avalanche, we wavelet-transformed the 
instantaneous configuration of interface heights, i.e. the 
set of {hi}, i = 1, . . . , L. As the avalanches progress, one 
gets a sampling of successive equilibrium configurations, 
which in turn provides us with an ensemble of the cor- 
responding wavelet coefficients. For each scale these are 
then translation-averaged, as explained above. 

In this case, the periodic boundary conditions imposed 
at the interface extremities are naturally consistent with 
those implicit in the wavelet transform, thus no potential 
mismatch arises. 

For comparison with BN simulation data, we gener- 
ated an artificial profile with ^ = 1.25, using the random 
midpoint displacement algorithm j3l|. Although earlier 
applications of wavelet transforms to fractional Brownian 
motion were restricted to < C < 1 in Ref.0, we found no 
technical impediments in going above that upper limit. 

It is known that profiles with C > 1 are rather 
smooth [H]. This is apparently at odds with the results 
to be expected from the force law, Eq. (HJ, from which 
the random locations of pinning centers would favor a 
rugged interface shape. Thus, it is worth looking at in- 
terface configurations in real space. One anticipates from 
Eq. idl) that the surface tension must play an important 
role in this context. Accordingly, we allowed k to vary 
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Figure 1: (Color online) Snapshots of typical interface con- 
figurations. All with 4096 sites, and periodic boundary con- 
ditions at the edges, (a)-(c): two-dimensional BN simula- 
tion, with varying surface tension (see Eq. ([ij), respectively: 
K = 1.0 (a), 3.0 (b), and 10.0 (c). (d): artificial profile with 
C = 1.25. 



by one order of magnitude. In Fig. [U one sees that on a 
fixed (system-wide) scale, the persistence trends charac- 
teristic of C > 1/2 are indeed reinforced by increasing k. 

One can have a quantitative understanding of the 
trends shown in Fig. [H with the help of wavelet trans- 
forms. The corresponding results are displayed in Fig- 
ure [21 where the horizontal axis is in units of inverse 
length scale, or "wavenumber" k = I /a. From scaling 
arguments [3|, the averaged wavelet coefficients W[/i](A:) 
for a self-similar profile are expected to vary as 



W[h]{k) 



^[(l/2)+Cl 



(8) 



A least-squares fit of a power-law dependence to the 
artificial-profile data for 64 < A: < 4096 gives C = 1.25(1). 
Such a central estimate and its uncertainty are both in 
line with corresponding results for < C < 1 0|. 

One sees that for BN data, ( ~ 1.25 holds only up to a 
crossover scale, which (as argued above) increases with k. 
This is illustrated in the inset of Fig. [JJ where a section, 
with 1/32 of the full length of the ragged k — 1.0 interface 
of Fig. [T] (a), is examined. On this scale, the profile is 
indeed much smoother than its parent. 

A fit of 64 < fc < 4096 data for k = 10.0 results in 
C — 1.19(3). This can be compared, e.g., with finite- 
size scaling estimates via Eq. ^ for the present model 
with K = 1.0, and a sequence of 400 < L < 1200 with 
O(IO^) configurations each, for which one quotes C = 
1.24(1) Is]. Eq. ^ reminds one that the latter method 
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Figure 2: (Color online) Double-logarithmic plot of averaged 
wavelet coefficients against wavenumber k. Symbols joined 
by continuous lines: wavelet transform of interface roughness 
data from two-dimensional BN simulation (interface dimen- 
sionality d' = 1). L = 4096, 10^ samples, with varying surface 
tension k (see Eq. (HJ). Crosses: wavelet transform of syn- 
thetic proffie with Hurst exponent — 1.25. L = 4096, 10"^ 
samples. Continuous line at bottom right has slope —1.75. 
Inset: section of length L' = 128 of typical profile for n = 1.0, 
illustrating interface smoothness on short scales (compare 
Fig. [11(a)). 



only considers fluctuations on short scales, thus in the 
present case it rightly captures the persistent behavior 
characteristic of that Hmit (at the expense of not being 
sensitive to the different trends that dominate the picture 
at larger scales). 

We conclude that the quantitative behavior exhibited 
by interface roughness in BN is likely to change when 
studied on varying length scales. Though a regime should 
exist, which displays close similarity to the Edwards- 
Wilkinson class of interface evolution problems, this 
should cross over to a more ragged picture on larger scales 
(the precise location of such change being determined by 
the interplay between quenched randomness and surface 
tension) . Wavelet transforms are thus a particularly suit- 
able method for the study of this problem, on account of 
the equal access that is provided to multiple length scales. 



IV. TIME SERIES AND POWER SPECTRA 



given in Eqs. H]) and 

Initially we consider the adiabatic limit of very slow 
driving. 

In experiment, the integrated signal J^^ V{t) dt is pro- 
portional to the magnetization change (number of up- 
turned spins) during the interval Ai. In the adiabatic 
approximation, a box-like shape is impHcitly assumed for 
each avalanche (i.e. details of the internal structure of 
each peak, as it develops in time, are ignored, on account 
of its duration being very short) , thus the instantaneous 
signal intensity (spike height) is proportional to the cor- 
responding avalanche size. 

As the signal is intermittent, there are signiflcant pe- 
riods (waiting times, WT) of no activity at all. Waiting 
time distributions for the adiabatic regime were exam- 
ined in Ref. Q. These were found to be rather flat, apart 
from (i) a sharp cutoff at the high end (related to the fl- 
nite cutoff in the avalanche size probability distribution) , 
and (ii) a number of peaks concentrated in a somewhat 
narrow region, which are associated to very frequent and 
small, spatially localized (i.e., non-critical) events involv- 
ing typically TV = 1 — 10 sites Q. 

We investigate the auto-correlations of two quantities, 
namely waiting times (WT), and avalanche sizes (i.e., 
BN spike voltages V) . For X = WT, V we calculate nor- 
malized, two-time connected correlations, averaged over 
t: 



{X{t)X{t + T)), 



- 1 



(9) 



For a system with L = 400, we have generated 2 x 10'' 
distinct time series of BN events. It is known [1] that, on 
account of the demagnetizing factor, size-size correlations 
are negative at short times, and decay with a characteris- 
tic relaxation time which (for this system size, and for the 
values of physical parameters used here) , is tq ~ 0.14 i. 
Thus, for each sample we calculated correlations in the 
range < t < R, R = 1.2, by scanning moving "win- 
dows" of width R along an interval of width lOi?. In 
preparation for ulterior wavelet analysis, the results were 
binned into N — 1024 equal-width bins. Our results are 
depicted in Fig. [H 

The exponential behavior of size data, noted ear- 
lier [1, is clearly discernible in the Figure even for 
T > 0.3, by which stage the signal-to-noise ratio has 
dipped to something close to unity. Waiting-time corre- 
lations initially seem to follow a similar exponential trend 
(with a time constant ~ 1/4 that for their size counter- 
part); however, a sharp "shoulder" develops at r « 0.1, 
signalling an abrupt end to the exponential regime. This 
indicates that negative waiting-time and size correlations 
have differing underlying causes. 



A. Introduction 

As explained above, owing to the assumed Hnear in- 
crease of applied field with time (in analogy with exper- 
imental setups), we shall express time in units of H as 



B. 



Waiting-time correlations 



Indeed, in calculating the correlations shown in Fig. [31 
the time separation r between any two waiting times is 




Figure 3: (Color online) Waiting-time (WT) and size corre- 
lations (see Eq. ([9])) against "time" in the adiabatic regime, 
for system with L = 400, 2 x 10* samples. Inset: absolute 
values of Gxij) on semi-logarithmic plot, same data range as 
in main Figure. 



considered to be the separation between their respec- 
tive starting moments (the same is done for size cor- 
relations, but it turns out to be of no further conse- 
quence, as avalanches are instantaneous in the adiabatic 
regime). This implies that the minimum separation be- 
tween two waiting times is the extent of the shortest of 
the two. Therefore, an effect arises at very short times 
T, which is the analogue of hard-core repulsion for sto- 
ichiometric problems in real space. Since the distribu- 
tion of waiting times is flat on a logarithmic scale [§| 
(thus P{WT) ^ 1/WT on a linear scale), and assuming 
waiting times to be uncorrelated (to be checked below), 
Eq. ^ gives |Gvkt(t)| « 1 — ar ~ exp(— ar) for r ^ 0. 

In order to eliminate this artifact, we then decided to 
index waiting times simply by their order of occurrence, 
thus (with j, k nonnegative integers) 



GwtU) = 



{WT{k)WT{k + j))k 
{WT{k))l 



- 1 



(10) 



In analogy with our earlier procedure, correlations were 
accumulated for j = 1, . . . , TV (iV = 1024) by generat- 
ing 20 independent series of 10 consecutive events; for 
each series we scanned moving "windows", each compris- 
ing iV -I- 1 events, i.e., TV waiting times. This time, the 
result was essentially flat noise, with no apparent short- 
time structure (see inset in Figure [4] below) . Therefore, 
further characterization must proceed via spectral analy- 
sis. We briefly recall how this can be done using wavelets. 

Assume one has 1//" noise. One calculates and 
wavelet-transforms the corresponding ensemble-averaged 



Figure 4: (Color online) Double-logarithmic plot of averaged 
wavelet coefficients against frequency /. Squares: wavelet 
transform of waiting-time autocorrelation data from two- 
dimensional BN simulation in the adiabatic regime, calcu- 
lated according to Eq. I|10p . L = 400, 20 independent series 
of 10 X 1024 waiting times. Crosses: wavelet transform of 
autocorrelations for synthetic 1//" noise, a = 1/2, L = 4096, 
5 X 10^ samples. A least-squares fit of 16 < / < 1024 data 
gives a = 0.51(1). Continuous line has slope —1/2. Inset: 
waiting-time correlations from BN simulation, calculated ac- 
cording to Eq. lflO|l . 



autocorrelations, and then translation-averages the re- 
sulting coefficients at each scale. Denoting the set of aver- 
aged wavelet coefficients by {VF[g](r)}, where {T} stands 
for the hierarchical set of wavelet timescales, and chang- 
ing the independent variable to "frequency" / = l/T, one 
expects from scaling [3. [2^: 



(11) 



For a > 1 this is derived immediately from Eq. |(8]), plus 
the exponent relation a = 1 -I- 2C (i,[2i|. Though for < 
a < 1 the scaling of cumulants of the noise distribution 
differs from that for a > 1, the basic scaling properties 
underlying Eq. (fTTI) remain valid 

Eq. (fTTj) can be tested with pure 1//" noise via the 
usual procedure of first producing a sequence of Gaussian 
white noise, Fourier-transforming that sequence, multi- 
plying the Fourier components by /~"^^ and then invert- 
ing the Fourier transform (23l.[33|. The resulting sequence 
is pure 1//" noise. An example with a = 1/2 is shown 
in Figure [4] below. 

Our results for BN are shown in Fig.jH Apart from the 
lowest frequency scale (which is not expected to fall in 
line with the rest, as it represents the most smoothed- 
out behavior [2]), the fiatness of the averaged coeffi- 
cients against varying scales strongly indicates that a = 



6 



(white noise), i.e., waiting times are indeed uncorrelated. 

The sequences of waiting-time correlation data of 
course need not be periodic. However, as seen above, 
they behave as random noise. Contrarily to, e.g., gen- 
erahzed brownian-motion profiles, such data are non- 
cumulative (i.e., they are not constrained in the fashion of 
consecutive positions of a random-walker, which cannot 
differ by more than one step length). Thus, the periodic 
boundary conditions implicit in the wavelet transform are 
not expected to introduce significant distortions in their 
analysis. 

C. Size correlations: adiabatic regime 

We now turn to the treatment of voltage data. For the 
adiabatic version of the interface model, of course only 
inter-avalanche voltage correlations can be evaluated. As 
mentioned above, the data in Fig.[3]are very well fitted by 
an exponential, with a "loading time" tq — 0.14(1). One 
then expects the Fourier power spectrum to be essentially 
fiat for / <C Tq"^, and to behave as 1/p for / ^ Tq^. 
This has indeed been found, e.g., in ref. 

The correlations to be wavelet-transformed are non- 
periodic and follow a clear baseline trend, therefore one 
needs to assess and eliminate potential distortions caused 
by: (i) using a periodic wavelet basis, and (ii) the baseline 
trend itself. 

In Fourier analysis, the standard way to deal with (i) 
is by zero-padding a region around the function to be 
transformed However, zero-padding does not work 
well when the function varies by orders of magnitude be- 
tween the extremes of the interval [2], as is the case here 
where only fiuctuations are left at the upper end. Tech- 
niques have been developed to remove the effects of peri- 
odic boundary conditions from wavelet transforms (i.e., 
to consider "wavelets on the interval") fsS]. These have 
very recently been translated into published computer 
code [3^, restricted to the Daub4 class. In the following, 
motivated especially by the need to address point (ii) , we 
propose a simplified approach based on detrending ideas. 
Combinations of wavelet decomposition and detrending 
have been investigated jsHi; however, the averaged coef- 
ficient analysis which is our main concern here has not 
been considered, except for some very simple cases (linear 
and quadratic drift |3|). 

We first illustrate how the averaged coefficients are af- 
fected by an overall exponential trend. Using the periodic 
Daub4 basis, we wavelet-transformed the size-correlation 
fitting function, Gyir) = -exp(-r/To). From Eq. 
one has: 

W[G^']{a,b) - — / V'a:6(a;)e--/^°dx . (12) 
Va J -00 

By changing variables, Eq. lfT2|) turns into: 

W[G^']{a,b) = V^e-'/"« f Vi;o(a;')e-""'/"Mx' . 

J —00 

(13) 



The first p = M/2 moments (starting at zeroth order) 
of Daubechies wavelets of order M vanish [2]. Thus, for 
M = 4 as is the case here, Taylor-expanding the expo- 
nential in the integrand of Eq. lfT3|) . one sees that the 
lowest-order non-zero term is proportional to a^/^, i.e., 

W[G^'](a,6) cxa5/2e-Vro _,.0(„7/2) ^^^^ 

We evaluated Gy* (r) at = 4096 equally-spaced points 
in the interval < t < 1.5, and wavelet-transformed it. 
For each hierarchical level j > 2, we plotted all 2^ wavelet 
coefficients, and found that they fall on the exponential- 
decay pattern of the original function, and (at the j-th 
hierarchical level) are proportional to 2"^-'/^, both fea- 
tures as predicted in Eq. |l4l) . except for the last two 
("wraparound" coefficients [2|). In order to fulfill the im- 
plicitly assumed periodicity of the original function, the 
latter coefficients take values ~ 10-' larger than the last 
preceding one (see an example for j = 5 in the inset of 
Fig. [U . Including these data in the coefficient-averaging 
procedure would introduce sizeable distortions (we did it, 
and found that the coefficients thus averaged behave as 
1//, which is in clear disagreement with the prediction 
of Eq. I|14p of a scaling power 5/2). 

To correct this artifact, we discarded the wraparound 
coefficients from the averaging procedure. Similar pro- 
cedures have been adopted elsewhere [sHi- As can be 
seen in Fig. \S[ this was enough to restore the expected 
behavior. Thus, point (i) above has been dealt with. 
We also wavelet-transformed Gy'(r) using the periodic 
Daubl2 basis. As expected, the coeficients behaved ap- 
proximately as a^'^/^ e"''/'^". The last four coefficients 
at each hierarchical level showed considerable increase 
against the exponential-decay pattern (as opposed to the 
last two for Daub4). In summary, as regards point (ii) 
we have shown that the most prominent feature of the 
wavelet transform (in the context of average wavelet co- 
efficient scaling), namely the Hurst-like exponent, of such 
a smooth function as the exponential fit is in fact basis- 
dependent. 

Thus, our simulational data must be detrended in or- 
der to eliminate distortions coming from the smooth 
baseline, which risk contaminating all scales. We did 
this by first subtracting the dominant exponential be- 
havior given by Gy*(T); for further refinement, we then 
removed some remaining non-monotonic mismatch via 
least-squares fit of a secondary adjusting function /(r) 
(a fourth-degree polynomial enveloped by a single ex- 
ponential), so G${t) = Gv{t) ~ G^*(t) - /(r). The 
result of wavelet-transforming the fully detrended cor- 
relations is depicted in Fig. [6l while the corresponding 
raw (detrended) data are shown in the inset of the same 
Figure (together with /(r), so one can have a quantita- 
tive estimate of how far the single-exponential fit goes 
to describe the un-detrended data). Note that /(r) has 
significant smooth variations on scales of 6t = 0.05 or 
longer, which translate into wavevectors A: < 32. We have 
wavelet-transformed partially-detrended data [i.e., with- 
out subtracting /(t)]. The respective averaged wavelet 
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Figure 5: (Color online) Main diagram: double-logarithmic 
plot of averaged wavelet coeflEicients against frequency /. 
Squares: wavelet transform of N = 4096 points of fitting 
function for size correlations, G|J'(r), for < r < 1.5. At 
each hierarchical level j > 2, the last two coefficients were 
omitted from the averages (see text). Continuous straight 
line has slope —5/2. Inset: semilogarithmic plot of (absolute 
value of) all 32 wavelet coefficients >V[Gv'](a, 6) [denoted by 
Wa(&)] against translation parameter b, at hierarchical level 

i = 5. 



coefficients are ~ 10 times larger than those for the fully- 
detrended curve for k < 16, and fall fast for increas- 
ing k: aX k = 64 the ratio is 1.4, and for A: > 64 both 
sets coincide to within less than 1%. So, failing to sub- 
tract /(r) introduces artificially large coefficients at large 
scales, which are not noise-related. 

Note that similar remarks apply here as in the earlier 
case of waiting-time correlations, namely, since Gy{t) is 
essentially noise around a horizontal basehne, the peri- 
odic boundary conditions implicit in the wavelet trans- 
form must not imply any significant distortion in our re- 
sults. 

The results exhibited in the main diagram of Fig. [6] 
strongly indicate that the detrended size correlations be- 
have as 1//" (white) noise. We defer discussion of this 
until the next Subsection, where departures from the adi- 
abatic regime are investigated. 



D. Size correlations: finite driving rate 

In order to discuss intra-avalanche correlations, one 
must introduce a finite driving rate 0, [U, so sep- 
arate events within the same avalanche can be ascribed 
to different instants in time. In fine with standard prac- 
tice [ll, [m, [13, m, [3§| our basic time unit is one lat- 



Figure 6: (Color online) Double-logarithmic plot of averaged 
wavelet coefficients against frequency /. Squares: wavelet 
transform of detrended size autocorrelation data, Gy{t) from 
two-dimensional BN simulation in the adiabatic regime. L = 
400, 2 X 10* samples. Inset: full lines: fully-detrended size 
correlations from BN simulation; dashed line: secondary ad- 
justing functon /(t) (see text). 



tice sweep, during which the external field is kept con- 
stant, and all spins on the interface are probed sequen- 
tially as described above. In the adiabatic regime, the 
external field is kept constant for the whole duration of 
an avalanche, i.e. for as many sweeps as it takes until 
no unstable sites are found along the interface. At finite 
driving rates, the field is increased by a fixed amount, 
henceforth denoted A, at the start of each sweep while 
an avalanche is taking place. Eventually, no more un- 
stable sites will be left, and then one proceeds as in the 
adiabatic regime, increasing the field by the minimum 
amount 6H necessary to start a new avalanche. In these 
"real" time units, the waiting time between the end of one 
event and the start of the next is then SH/ A; however, in 
order to produce meaningful comparisons, especially be- 
tween data acquired in the adiabatic and non-adiabatic 
regimes, it will be useful to keep referring to the "ab- 
solute" scale given by the applied field H itself, which 
unequivocally locates events along the hysteresis cycle. 

As A grows, the intermittent character of events is 
gradually lost as more and more avalanches coalesce 
and one eventually crosses over to a regime in which the 
interface is fully depinned, i.e. it moves at non-zero av- 
erage speed. 

In Fig. [7] we show autocorrelations for driving rates still 
within the intermittent regime, compared with those for 
the adiabatic Hmit. The most significant change upon 
increasing A is the effective loss of negative short-time 
correlations. In fact, this represents an excess of posi- 



8 




Figure 7: (Color online) Normalized two-time correlations 
(averaged over t) {V{t)V{t + T))/{V{t)Y - 1 from two- 
dimensional BN simulation, for system with L — 400, and 
driving rates A as given in key to symbols (A = corre- 
sponds to adiabatic limit). "Time" is given in applied field 
units, i.e. "absolute" scale (see text). 



tive intra- avalanche contributions, on top of the negative 
inter-avalanche terms (and some intra-avalanche ones as 
well) which still exist for non-zero A (on account of 
the demagnetizing factor). Positive reinforcements arise 
mostly because, when many sites are overturned during 
one lattice sweep, that same number of new sites will be 
probed by the interface. For each new site, the quenched 
randomness term in Eq.[T]may, or may not, contribute to 
further motion with roughly equal chances. By contrast, 
at a site which remains pinned during one sweep, the in- 
terface stands fewer chances of getting unstuck, as the 
contribution from the randomness term is kept constant; 
depinning of such a site is more likely to happen if the 
field is substantially increased, i.e., during a subsequent 
avalanche. 

We detrended the A 7^ data of Fig. [7] by similar pro- 
cedures to those used earlier for A = 0. The main differ- 
ence was that detrending was done in a single stage, fit- 
ting /(r) described in Subsec. ITV CI to the raw data, and 
then subtracting the least-squares-fit from the original 
data. The results of wavelet-transforming the detrended 
data are shown in Figure [8l 

One can see that, as opposed to the adiabatic regime, 
data for finite driving rates clearly exhibit a downward 
trend for a range of intermediate frequencies, spanning 
3-4 hierarchical levels, and which is characterized by an 
approximate 1//^'^ behavior (the straight line segments 
in the Figure have slope —1.5). Furthermore, with the 
"absolute" frequency / given in inverse applied field units. 



Figure 8: (Color online) Double-logarithmic plot of averaged 
wavelet coefficients against frequency /, from wavelet trans- 
form of detrended size autocorrelation data, Gy{t), for two- 
dimensional BN simulations of system with L = 400, and 
assorted driving rates A. Key to symbols is the same as in 
Fig. [7] (A = corresponds to adiabatic limit). Frequency 
is given in inverse applied field units, i.e. "absolute" scale 
(see text). Plots successively shifted downward by a factor 
of 10 on vertical scale, to avoid superposition. Straight line 
segments mark subsets of A 7^ regime where approximate 
1//^-^ behavior holds. 



and A given in units of applied field change per unit time, 
dimensional arguments show that /' = / A is the "nat- 
ural" frequency variable (i.e., inverse "real" time). This 
is shown more clearly on a scaling plot. Figure [9l where 
use of /' as the independent variable causes the 1//^'^ 
sections of all A data to collapse. 

Given that, in these slow- (but non- adiabatic) driv- 
ing regimes, avalanche coalescence comprises only a small 
fraction of events Q, one can say that approximately the 
same sequence of avalanches occurs for all A investigated 
here, only at different "real" paces. Since the "real" time 
interval between consecutive avalanches is Sh/ A, and as- 
suming Sh to be the same, for different values of A, be- 
tween two given avalanches (for the reasons just men- 
tioned), one sees that inter-avalanche correlations will 
shift to higher "real" frequencies as A grows. On the 
other hand, within a given avalanche, two sub-events 
separated by a given number of lattice sweeps are (by 
definition used in the simulation) separated by the same 
"real" time interval, thus their correlations are not shifted 
in "real" frequency for varying A. Therefore we conclude 
that the collapsing sections of the scaling plot correspond 
mainly to intra-avalanche correlations. 

First-order (Haar) spectra of experimental BN data 
show that, for Fe2iCo64Bi5 the high-frequency section 
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Figure 9: (Color online) Double-logarithmic scaling plot of 
averaged wavelet coeflEicients against "natural" frequency /' = 
/ A, from wavelet transform of detrended size autocorrelation 
data, Gv{t), for two-dimensional BN simulations of system 
with L = 400, and assorted driving rates A 7^ 0. Key to 
symbols is the same as in Figs. [7] and [8] Full straight line has 
slope —1.5 . 

falls initially as /~^'^, and then crosses over to f~^'^, 
while for FeSi the decay is with /"^-^^ Though 
the exponent values in both cases are not too dissimilar to 
the one found here, analysis of higher-order spectra 
leads to a more nuanced picture. For Fe2iCo64Bi5, it 
is found that most of the power in the high-frequency 
range comes from intra-pulse correlations [l9| , similar to 
our conclusion above, whereas for FeSi the conclusion 
was that the high-frequency power is mainly connected 
to the inter-pulse sort [131 • Therefore it would appear 
that the dynamics of the present model is closer to that 
of BN in materials like Fe2iCo64Bi5 than in FeSi. 

V. DISCUSSION AND CONCLUSIONS 

We have discussed the application of wavelet trans- 
forms to the description of both real-space and time-like 
properties of an interface model, which is used for the de- 
scription of Barkhausen noise in soft ferromagnets. Most 
of our calculations involved the scaling properties of posi- 
tional averages of wavelet coefficients, taken at each hier- 
archical (size) level, as first proposed in Ref. S In some 
instances we showed that direct analysis of individual 
coefficients was called for, in order to unravel artificial 
effects which would otherwise distort our aggregate re- 
sults. 

Here we considered the d = 2 version of the model 
(thus the interface dimensionality is d' = 1), mainly in 



the adiabatic limit of very slow driving, for which the 
sudden "avalanches" of domain wall motions are consid- 
ered to occur instantaneously. In Subsection lIVDi we 
extended our study to finite driving rates, in order to 
analyze intra-avalanche correlations 

Our investigation of real-space aspects consisted in the 
evaluation of the characteristic interface roughness expo- 
nent C.. On scales shorter than a crossover length (which 
turns larger as the intensity of surface tension grows), 
we get C = 1-20(3), close to C = 1-24(1), derived by 
other methods for the same model Q, and also to as- 
sorted estimates for quenched Edwards- Wilkinson sys- 
tems |23,[2i,S[33|, which give C ^ 1-25. 

Turning to time series, in Subsection IIV Bl we showed 
that a proper indexation of the sequence of waiting-times 
between avalanches is crucial, in order to avoid artificial 
short-time negative correlations. Procedures similar to 
that used here, namely, indexing waiting times simply 
by their order of occurrence (instead of using the start- 
ing time of each interval) , have been used consistently in 
the context of self-organized criticality scaling Our 
final result (see Fig.|4]) was that the correlations between 
waiting times are white noise, i.e., these quantities are 
fully uncorrelated. Going back to the rules of interface 
motion, and to Eqs. |(T]) and ([2]), one sees that this is a sig- 
nature of the quenched-randomness term u{xi,hi). This 
fact is in contrast to the behavior of size correlations, 
which are strongly infiuenced by demagnetization [4, 9]. 

In Subsection IIV C\ we started from the known fact 
that, in the adiabatic regime, size-size correlations are 
negative at short times, and decay approximately as an 
exponential 0, 0|. By direct analysis of (non-averaged) 
wavelet coefficients, we illustrated practical ways to deal 
with artifacts introduced by the periodicity of the wavelet 
basis used. It turned out that the smooth baseline func- 
tion, to which noise data are fitted, can introduce distor- 
tions at all levels of the wavelet transform. Furthermore, 
such distortions are non-universal in the sense that they 
depend on the wavelet basis. Thus, in order to obtain 
meaningful results from averaged wavelet coefficients, one 
must fully detrend the raw data. Once we did so, we 
found strong indications that the detrended size corre- 
lations behave as white noise (see Fig. [6]). This is ap- 
parently at odds with earlier (Fourier) power-spectrum 
results (see, e.g., Ref. @ and references therein), which 
would lead one to expect 1/p behavior, at least for high 
frequencies. However, the derivation of the latter result 
(e.g., by direct integration) fully takes into account the 
exponential basehne shape, thus one is referring to a dif- 
ferent object. Here, as explained above, we are deafing 
with detrended data. 

Finally, in Subsection IIV T)\ we considered size corre- 
lations against time in non-adiabatic regimes (but well 
within the driving-rate range where intermittency still 
holds 0). For driving rates A spanning one order of mag- 
nitude, we found rather well-defined frequency intervals 
for which detrended correlations behave as /~", a « 1.5, 
By changing variables from "absolute" to "natural" fre- 
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quency, we found that said intervals collapse together, 
which indicates that they pertain to intra-avalanche cor- 
relations. Rather than attaching much significance to 
the numerical value of the power-law exponent (since 
the shortness of the interval along which such behavior 
holds prevents one from doing so), one must emphasize 
the good degree of curve collapse exactly in that sec- 
tion, and only there. This indicates that this section is 
the "special" one, i.e. it corresponds to the frequency 
range along which universal (driving-rate independent) 
properties hold. Furthermore, our considerations lead- 
ing to the conclusion that such scaling behavior reflects 
intra-avalanche correlations are completely independent 



of the analysis of higher-order power spectra experimen- 
tal data, carried out in Refs. O, [13, and which leads to 
the very same conclusion as regards BN in samples of 
Fe2iCo64Bi5. 
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